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We investigate the relationship between synchronized muscle activity and tremor time series in (en- 
hanced) physiological tremor by cross-spectral analysis. Special attention is directed to the phase 
spectrum and its possibilities to clarify the contribution of reflex mechanisms to physiological tremor. 
The phase spectra are investigated under the assumptions that the EMG synchronization was caused 
by a reflex, respectively a central oscillator. In comparison of these results to phase spectra of mea- 
sured data we found a significant contribution of reflexes. But reflexes only modify existing peaks in 
the power spectrum. The main agents of physiological tremor are an effterent pace and the resonant 
behavior of the biomechanical system. 



1 Introduction 

The contribution of reflex loops to (enhanced) physiologi- 
cal tremor is still a matter of debate (Elble and KoUer 
1990). Most experiments investigating this topic used 
some kind of manipulation of the system. For example, 
Lippold (1971) cooled the limb in order to reduce the ve- 
locity of nerve conduction. Rack et al. (1978) examined 
the reflex response to enforced sinusoidal movement of the 
limb and Young and Hagbarth (1980) investigated several 
manoeuvres that affect the stretch reflex to enhance phys- 
iological tremor (ETP). We intend to investigate the role 
of reflex mechanisms by cross-spectral analysis of surface 
electromyogram (EMG) and accelerometer data (ACC) 
without perturbing the system externally. Former appli- 
cations of this technique to physiological tremor are given 
in (Fox and Randall 1970; Pashda and Stein 1973; Elble 
and Randall 1976; laizzo and Pozos 1992). These authors 
interpreted the phase spectrum at single, fixed frequen- 
cies in terms of delays between the EMG and the AC.C 
We examine the whole phase spectrum and investigate 
whether the phase spectrum can help to clarify the role 
of reflexes in EPT. 

Typical EPT usually shows two significant peaks in the 
ACC power spectrum. These peaks exhibit a different be- 
havior if the hand is loaded. While often indistinguishable 
located around 10 to 12 Hz if the hand is not loaded, one 
peak moves to lower frequencies under increasing loads. 
This peak is called mechanical peak since it depends on 
the mechanical properties of the musculosceletal system. 
An EMG peak at the same frequency shows up occasion- 
ally. The frequency of the second peak does not vary 
under moderate loading and is always accompanied by a 
peak in the EMG power spectrum (Homberg et al. 1987; 
Elble and KoUer 1990). This peak is called neurogenic. 



It is still controversial whether these EMG synchroniza- 
tions result from reflex loops, by chance synchronization 
or from a central oscillator (see Allum et al. 1978; Stiles 
1980; Christakos 1982; Allum 1984; Elble and KoUer 1990 
for a review). The two latter hypotheses can not be dis- 
tinguished by cross-spectral analysis. 

For both hypotheses of reflex mechanisms and non- 
reflex mechanisms, we calculate the power spectrum, the 
cross-correlation function and the phase spectrum. In the 
case of the non-reflex hypothesis, it is possible to derive 
these quantities analytically. In the case of reflexes, some 
kind of nonlinearity has to be introduced into the feedback 
system which hinders an analytical treatment. Based on 
a stochastic feedback system for the effects of reflexes we 
perform simulation studies to understand the behavior of 
the cross-correlation function and the spectra. This model 
is a generalization of the deterministic model introduced 
by Stein and OguztoreH (Stein and OguztoreH 1975, 1978; 
OguztoreU and Stein 1975, 1976, 1979; Bawa et al. 1976a, 
1976b). 

It should be mentioned that the well developed theory 
of Volterra - Wiener kernel estimators, which has been 
successfully applied to model physiological feedback sys- 
tems (Marmarelis 1989) can not be applied to the problem 
under investigation since it requires the input data of the 
system, which are not available here. 

We do not intend to derive a quantitative model of the 
reflex mechanism. We will use the phase spectra as an 
indicator for the presence of reflex contributions by com- 
paring the predicted phase spectra under both hypotheses 
with those estimated from measured data. For both cases 
the mechanical properties of the hand are modeled by 
an autoregressive process of order to two (AR[2]) covered 
by observational noise (Stiles and Randall 1967; Randall 
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1973; Rietz and Stiles 1974; Elble and Roller 1990; Gan- 
tert et al. 1992). The synchronized (possible) non-reflex 
part of the rectified EMG is also described by an AR[2] 
process. We claim no physiological significance for this 
EMG model; it is used solely as a description of the data 
and an approximation to the true statistical behavior of 
this part of the EMG. Since we only investigate second 
order statistical properties of the signals this description 
is adequate if the power spectra of the model and the 
data are similar. This model is justified by the over-all 
consistent description of the phenomena observed in the 
data. 

In a companion paper (Timmer et al. 1998), we have de- 
scribed the recording conditions, the mathematical meth- 
ods and their application to EMG and ACC data for cases 
where the EMG does not show significant synchroniza- 
tion. 

This paper is organized as follows. In the next sec- 
tion we discuss the characteristic features of the phase 
spectrum between EMG and ACC assuming that the syn- 
chronized EMG does not result from a reflex mechanism. 
In Section ^ we give the corresponding results assuming 
a significant contribution of reflex loops. In Section ^ we 
give the results for 35 data sets recorded from 19 subjects. 

2 Non reflex mechanisms 

In this section, we show the expected behavior of the dif- 
ferent kinds of spectra assuming that the synchronized 
EMG activity is not caused by any reflex loop. Such a 
synchronization can result from motor units flring asyn- 
chronously at similar rates (Christakos 1982) or some cen- 
tral oscillator. Cross-spectral analysis will not be able 
to distinguish between peripheral and central non-reflex 
mechanisms. We first show results of simulations and then 
applications to measured data. 

2.1 Simulations 

If the synchronization of the EMG is not caused by a re- 
flex, the position of the hand, x{t) (the ACC is obtained 
by differentiating x{t) twice), is described by an AR pro- 
cess driven by colored noise. Denoting uncorrelated EMG 
activity by i^(t), the synchronized EMG by y{t), and a 
possible time delay of the effect of the EMG on the ACC 
by At, the model reads: 

y{t) = b,y{t-l) + b2y{t-2) + e{t) (1) 
x{t) = aix{t - 1) + a2x{t - 2) (2) 
+iy{t - At) + y{t - At) . 

The parameters ai and a2 describe the physical proper- 
ties of the hand/muscles system, i.e. resonance frequency 
l/T and damping rate r. hi and 62 were chosen in or- 
der to simulate the desired power spectra of the EMG. 
Observational noise was added to obtain the measured 



EMG and ACC. The power spectrum of the ACC is given 
by (15) of the companion paper (Timmer et al. 1998). 
The phase spectrum does not depend on the properties 
of the EMG activity. Fig. | gives the results for the case 
that the EMG peak shows a higher frequency (14 Hz, 
hi = 1.8924,62 — —0.9769) than the resonance frequency 
of 10 Hz [ai = 1.9401, 03 = -0.9835) of the hand. This 
situation is similar to that often found in physiological 
tremor when the hand is loaded. Fig. |^a shows the power 
spectra of EMG and ACC, Fig. |l|b the coherency spec- 
trum and Fig. ^ the phase spectrum which is displayed 
27r periodically for a range of ±37r. Fig. |l|d the cross- 
correlation function. In Fig. ^, results for the opposite 
peak frequency location (ai = 1.8318, 02 — —0.9704, hi = 
1.8924,62 = -0.9769 ) are displayed. This situation is 
similar to the finger tremor with a resonance frequency of 
about 20 Hz, whereas the EMG spectrum remains usually 
unchanged. Note, that the phase spectrum only depends 
on the resonant system. 

Phase spectra like those in Fig. |l| und ^ have been inter- 
preted in earlier publications only at the frequency wq of 
maximum coherency in terms of a delay At (Pashda and 
Stein 1973; Elble and RandaU 1976) by At = $(u;o)/wo- 
Obviously, such an interpretation might be misleading, 
since no delay is included in the process. The frequency 
of maximum coherency is determined by the frequency 
dependent signal to noise ratio ((27) of the companion 
paper) und usually occurs at the peak frequency of the 
EMG. The phase spectrum evaluated at the peak fre- 
quency of the EMG only contains the information whether 
the resonance frequency of the hand/muscle system is lo- 
cated above or below the EMG peak frequency. 

The cross-correlation function is difficult to interpret. 
The period of the cross-correlation function is mainly de- 
termined by the frequency of highest coherency and mod- 
ulated by frequencies of medium coherency. Obviously, 
the time-lag where the cross-correlation function shows 
its maximum may, in general, not be interpreted as a time 
delay between the processes. The magnitude of the cross- 
correlation function is similar for positive and negative 
lags. Especially for negative time lags, it exceeds signif- 
icantly the 5 % significance level of ±1.96A^~^ of zero 
cross-correlation derived under the assumption that at 
least one of the independent processes is white noise. The 
reason is that the autocorrelation function of both pro- 
cesses and the true cross-correlation function enters the 
correlation structure of the estimated cross-correlation 
function here. Thus, large values of the cross-correlation 
function for negative time lags, suggesting an effect from 
the ACC on the EMG, do not necessarily indicate the 
involvement of a reflex loop. 

2.2 Application to measured data and 
modeling 

Fig. H shows the results of an analysis of a physiological 
hand tremor showing an EMG synchronization at approx- 
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Figure 1: Results of a simulations study, a: power spectra 
(EMG: dashed line, ACC: solid line), b: coherency spectrum, 
the straight line in represents the 5% significance level for the 
hypothesis of zero coherency, c: phase spectrum with 95% con- 
fidence interval, plotted 2n periodically, d: cross-correlation 
function with 5% significance levels for the hypothesis of zero 
cross-correlation assuming that at least one of the processes is 
white noise. The peak frequency of the EMG is higher than 
the resonance frequency of the hand. The step in the phase 
spectrum is determined by the latter. 



imately 10.5 Hz. It resembles to Fig. |^. The behavior of 
the coherency can be explained as an effect of observa- 
tional noise by (27) of the companion paper. In order to 



Figure 2: Results of a simulations study, a: power spectra 
(EMG: dashed line, ACC: solid line), b: coherency spectrum, 
c: phase spectrum, d: cross-correlation function. The peak 
frequency of the EMG ts lower than the resonance frequency 
of the hand. The step in the phase spectrum is determined by 
the latter. 



show that these results can be interpreted in the frame- 
work of linear models we fitted a model of seven param- 
eters according to (|l|,^ to the empirical data. The seven 
parameters are periods and relaxation times describing 
the EMG and the mechanical properties of the hand, 
the variance of the observational noise covering EMG 
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and ACC and the amount of additional unsynchronized 
EMG activity. By Equations (2,3,4,5,13,15,22,23,27) of 
the companion paper these parameters can be estimated 
from the data. Fig. ^ shows the results of a simulation 
of the fitted model. All quantitative details of Fig. ^ are 
reproduced by the model. 
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Ftgure 3: Physiological tremor with EMG synchronization, 
a: power spectra (EMG: dashed line, ACC: solid line), b: 
coherency spectrum, c: phase spectrum, d: cross-correlation 
function. 
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Figure 4: Realization of the model fitted to the data of Fig. 3. 

a: power spectra (EMG: dashed line, ACC: solid line), b: 

coherency spectrum, c: phase spectrum, d: cross-correlation 
function. 



3 Reflex mechanisms 

Pioneered by the work of Lippold (Lippold 1957, 1970) 
various experiments were performed to clarify the effect 
of reflex mechanisms on physiological tremor. In most of 
these experiments a mechanical or electrical stimulus was 
applied to measure the response of the muscle (Rack et 
al. 1978, Hagbarth and Young 1979, Matthew 1994). By 
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cross-spectral analysis we studied the role of reflexes in 
physiological tremor without perturbing the system, i.e. 
the stimulus for the reflexes is the tremor itself. 

3.1 Simulations 

To investigate the phenomena that might be induced by 
reflex mechanisms, we use a stochastic feedback system 
which applies a sigmoidal nonlinearity describing the ac- 
tivation function of the motoneurons. The effects of such 
an activation function in a deterministic framework sim- 
ulating spinal as well as cortical feedback delays were in- 
troduced by Stein and Oguztoreli (Stein and Oguztoreli 
1975, 1978; Oguztoreh and Stein 1975, 1976, 1979). The 
reflex model of Stein and coworkers includes the muscle 
length and velocity. We use the tremor amplitude and 
velocity. Because of the small amplitude these quantities 
are proportional. In order to study the influences of reflex 
mechanisms on the spectra, we used the following model. 
Denoting the reflex induced EMG by y{t), the unsynchro- 
nized EMG activity by i'{t), the movement of the hand 
by x{t), a possible time delay of the effect of the EMG 
on the hand by At, and the reflex loop delay by 6t we 
receive: 

y{t) = Kf{x{t-St),x{t-dt)) (3) 
x{t) = aix{t - 1) + a2x{t - 2) (4) 
+iy{t - At) + y{t - At) . 

For the function /(.) wc chose the tangens hyperbolicus 
with gain n; oi and 02 describe the physical properties of 
the musculosceletal system. x{t — St) models the length 
dependent and i{t—6t) the velocity dependent reflex. Ob- 
servational noise was added to both processes. Its stan- 
dard deviation was 10 % of the standard deviation of the 
noise free data. Note, that form a nonlinear stochas- 
tic delay differential equation. 

Due to the nonlinear activation function, in general, the 
spectra can not be calculated analytically. We performed 
simulation studies in order to investigate the behavior of 
the spectra. The unique result for various settings of the 
parameters is, that 

• the EMG power spectrum shows a peak at the same 
frequency as the ACC power spectrum. The former 
might be hidden by the observational noise. 

• in accordance with Oguztoreli and Stein (1975) the 
frequency of the oscillation is a complex function of 
the delays and the properties of the hand and mus- 
cles. Fig. ^ displays the reflex induced frequency 
shift in dependence on the delay of the reflex loop. 
Fig. ||b shows this shift as a function of the charac- 
teristic period of the mechanical system. 

• the phase spectrum is given by the mechanical prop- 
erties of the hand/muscle system represented by ai 
and 02, i.e. not dependent on the reflexes. 



• the peaks of the ACC are sharper in presence of a 
reflex than without, in accordance with experimental 
findings (Marsden 1984). 
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Figure 5: Change of peak frequency due to the reflex, a: De- 
pendence on the delay of the reflex loop for a fixed character- 
istic period (100 ms) of the mechanical system, b: Dependence 
on the period of the mechanical system for a fixed delay of the 
reflex loop of 30 ms. 

In general a non-rcflcx and a reflex component might 
show up in the EMG. For this case, the models (|l|,||) and 
(||,||) have to be combined to: 

V{t) = Kf[x{t-5t),x{t-5t)) (5) 
+hiy{t-l) + b2y{t-2) + t{t) 

x{t) = aix{t -l) + a2x{t - 2) (6) 
-f //(t - At) + y{t - At) . 

A representative example for this model without any 
reflex loop (k = 0), i.e. model (P|), is given in Fig. |6^-c. 
The non-reflex EMG peak at 10 Hz {bi = 1.9401,62 = 
—0.9835) causes an ACC peak at the same frequency. 

The ACC peak located at 5 Hz refers to the resonant 
behavior of the hand (oi = 1.9780,02 = -0.9889). As 
mentioned in Section 3.3 of the companion paper (Timmer 
et al. 1998) the phase spectrum is determined by the 
resonance properties of the hand/muscle system. Fig. 
f displays the results for a delay of the reflex loop of 35 
ms corresponding to a possible segmental stretch reflex 
(Noth et al. 1985). k was chosen equal to minus one. 
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Figure 6: Results of a simulation study, a: EMG (dashed) and ACC (solid) power spectra in the case of no reflex, 
corresponding phase spectrum c: coherency spectrum, d: EMG and ACC power spectra in the case of a segmental reflex, 
corresponding phase spectrum, f: coherency spectrum. 



The ACC peak caused by the resonant behavior of the 
hand is shifted to a higher frequency of 6.3 Hz. But 
the phase spectrum is invariant and still represents the 
characteristic features of the mechanic properties of the 
hand/muscle system that have shown up in the ACC 
power spectrum of Fig. |6^. Therefore, a significant differ- 
ence in the peak frequency estimated from the phase spec- 
trum and the power spectrum indicates a contribution of 
the reflex to the oscillation. There might be parameter 
constellations as shown in Fig. || where a reflex does not 
lead to any shift of the resonance frequency. Thus, we 
modify in the empirical study the mechanical properties 
of the system by loading the hand with different weights 
to increase the probability of detecting a possible contri- 
bution of a reflex loop. 

The peak frequency is estimated from the phase spec- 
trum by fitting the expected theoretical phase spectrum 
to the empirical one as described in Section 3.3 of the 
companion paper. The model contains a possible time 
delay At, the parameters oi and 02 of the AR[2] process 
describing the properties of the hand/muscle system. Fur- 
thermore, the twofold differentiation to obtain the ACC 
from the movement x(t) of the hand is considered. There- 
fore, the theoretical phase spectrum consists of the sum of 
(20,22) of the companion paper and an offset correspond- 
ing to (21). The peak frequency is determined by (2,3,14) 
of the companion paper from the estimated parameters 
oi and a2- 

Note that the coherency at the frequency of the reflex 
induced EMG peak increases due to the larger signal to 
noise ratio in accordance with (27) of the companion pa- 
per. 

3.2 Application to measured data 

Fig. and Fig. || display 3 s segments of data recorded 
from two subjects with FTP. Fig. ^ and |l^ show the re- 
sults for these two examples. In Fig. ^a the mechanical 
peak of the ACC at 6.2 Hz is accompanied by a peak in 
the EMG, located at the same frequency. The neurogenic 
EMG peak around 13 Hz causes a shoulder in the ACC 
power spectrum according to (15) of the companion pa- 
per. Because of the increased signal to noise ratio at the 
peaks, the coherency in Fig. ^b shows significant values 
there. The phase spectrum in Fig. ^c exhibits the ex- 
pected shape in dependence on the resonance properties 
of the hand/muscle system. The frequency of the mechan- 
ical system estimated by fitting the analytical phase spec- 
trum (20,21,22) of the companion paper to the estimated 
phase spectrum is 7.4 ± 0.1 Hz. The peak in the ACC 
power spectrum, representing the mechanical resonance 
frequency, is located at 6.2 ± 0.05 Hz (the calculation of 
the confidence regions is explained in Section 3.3 of the 
companion paper). Therefore, as the probability that the 
two frequencies are the same is less than 0.01, we conclude 
that the reflexes contribute signiflcantly to this tremor. 
In Fig. Il^a the mechanical peak at 6.9 Hz is not accom- 




Figure 7: Enhanced physiological hand tremor, a: Accelera- 
tion of the hand, b: rectified EMG 
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Figure 8: Enhanced physiological hand tremor, a: Accelera- 
tion of the hand, b: rectified EMG. 



panied by a significant peak in the EMG. The neurogenic 
peak at 11 Hz causes a further peak in the ACC power 
spectrum. The frequency of the mechanical system es- 
timated from the phase spectrum is 7.7 ±0.1 Hz. An 
estimation in the ACC spectrum yields to 6.9 ± 0.1 Hz. 
Thus, again, a refiex mechanism is involved (p < 0.01). 
Due to the reflex contribution one would expect a peak in 
the EMG spectrum located at the mechanical ACC peak. 
It might be not observable since the induced EMG activ- 
ity does not exceed the observational noise significantly. 

Because of the nonlinear structure of model (||J^) it can 
not be fitted directly to the data. 
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Figure 9: Analysis of data shown in Fig. 7. a: power spectra 
(EMG: dashed line, ACC: solid line), h: coherency, c: phase 
spectrum. The mechanical peak frequency estimated from the 
phase spectrum is 7.4 ± 0.2 Hz is not consistent with that esti- 
mated from the ACC power spectrum of 6.2 ±0.1 Hz. 



Figure 10: Analysis of data shown in Fig. 8. a: power spectra 
(EMG: dashed line, ACC: solid line), b: coherency, c: phase 
spectrum. The mechanical peak frequency estimated from the 
phase spectrum is 7.7 ± 0.2 Hz is not consistent with that esti- 
mated from the ACC power spectrum of 6.9 ± 0.2 Hz. 



4 Results 

In this section we report the results obtained by apply- 
ing methods introduced to 57 recordings of physiological 
tremor and EMG measured from 19 subjects. For each 
subject the data were recorded without loading the hand, 
with a load of 500 gr and 1000 gr in order to modify the 
resonance frequency of the hand. 

As described in Section 3.1 above and Section 3.3 of 
the companion paper, we fitted theoretical phase spectra 
to the empirical phase spectra to estimate the frequency 
of the mechanical system. If no reflex mechanism was 
involved this frequency would be located at the peak fre- 
quency in the ACC power spectrum. 

To obtain a good fit of the phase spectra it was always 
necessary to correct for the spurious time delay which is 



present because we modeled a continuous time processes 
by a discrete time model due to (23) of the companion 
paper. No further time delay was found in the data. In 
24 of the 57 recordings the errors in the phase spectrum 
were too large to obtain reliable results. This is caused 
by a poor coherency due to a small signal-to-noise ratio. 
Thus, only 35 data sets were included in the analysis. 

The differences between the frequencies estimated from 
the phase and the ACC power spectrum for different load- 
ings of the hand are given in Fig. |l^. The error-bars rep- 
resent the 2(7 confidence interval. Under the hypothesis 
that the reflex does not contribute to the process only two 
of the data points should not be consistent with a zero dif- 
ference at a signiflcance level of 5 %. 19 data points are 
signiflcantly different from zero. On the one hand this 
result show that reflexes play a role in ETP. On the other 
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Figure 11: Differences in peak frequencies estimated from 
the phase spectrum and the ACC power spectrum for different 
loadings (O: no loading, □; 500 gr., A: 1000 gr.) for 19 
subjects with EPT. 

hand it shows that a simple segmental stretch reflex is not 
suflicient to explain this finding. This is indicated by the 
behavior of the peak frequency differences in dependence 
of the loading. Simulation studies like those presented in 
Fig. ^ show that the difference should become more posi- 
tive under loading. The results allow no clear cut decision 
on this but the different trend suggests the involvement 
of more complex reflex structures in the process, maybe 
a combination of different reflex loops. 

5 Summary 

We investigated the relation between synchronized EMG 
and ACC of physiological tremor by cross-spectral analy- 
sis. The phase spectrum between EMG and ACC depends 
only on the mechanical properties of the hand/muscle sys- 
tem, i.e. the driven part of the system, but not on the 
characteristics of the driving force. The behavior of the 
coherency with respect to the peaks of EMG and ACC 
power spectra can be explained by its dependence on the 
observational noise. Usually the maximum coherency is 
located at the frequency of the EMG peak. Therefore, 
the value of the phase spectrum at the frequency of the 
EMG peak may not be interpreted as a delay. It only 
provides information about the relative location of the 
mechanical resonance frequency and the synchronization 
frequency in the EMG. The time delay of reflex loops can 
shift the resonance peak frequency, but do not affect the 
phase spectrum. Thus, in both cases with and without a 
contribution of reflexes, the phase spectrum provides the 
information about the frequency of the mechanical peak 
that would show up in the absence of any reflex. There- 
fore, a signiflcant difference in the peak frequency of the 
ACC power spectrum and that calculated from the phase 
spectrum gives evidence for a contribution of reflexes to 
the tremor. In 35 time series recorded from 19 subjects 



with ETP, we found clearly that reflexes contribute to 
this tremor. The sign of the difference was found to be 
opposite to that expected for segmental reflexes from sim- 
ulation studies suggesting that more complex structures 
are involved. 

There is no evidence in the data that reflex loops pri- 
marily cause the tremor. They alter the frequency, relax- 
ation time and amplitude of existing oscillations to some 
degree. This holds for the mechanical peak as well as 
for the so called neurogenic, i.e. centrally evoked, peak. 
Therefore, the primary cause of ETP is the resonant be- 
havior of the hand and a synchronized EMG activity that 
is either generated centrally or due to the recruitment 
strategy of motoneurons. 
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